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A theoretical analysis is presented of the relaxation processes which lead to the magnetization 
of gaseous hydrogen atoms in a strong magnetic field. Experiments are described which demonstrate 
that the spin-lattice relaxation time in a pure H 2 carrier can be of the order of 10 ms. From the effects 
of traces of 2 on the hydrogen atom EPR spectrum, a H-0 2 spin-exchange collision cross section of 
(8.0 ±2.5) X 10- ,(i cm 2 was derived. 

Keywords: Electron spin resonance; hydrogen atoms: spin-exchange: spin-relaxation. 



1. Introduction 

It has been demonstrated that gaseous nitrogen 
atoms can have a spin-lattice relaxation time 7\ of 
the order of several milliseconds [1J. 1 In high speed 
discharge-flow systems with EPR atom detection, this 
T\ can be longer than the transit time of the atoms 
between the discharge and the EPR microwave cavity. 
Anomalously small EPR absorption intensities will 
arise in this situation because the atoms entering the 
cavity will not be completely magnetized. Decreasing 
the gas flow rate or decreasing T\ by adding a para- 
magnetic impurity like 2 will increase the magnetiza- 
tion at the cavity. The EPR signal will thus increase 
even if the atom concentration in the cavity remains 
constant. To measure atom concentrations quantita- 
tively by EPR under such conditions one must know 
how far the magnetization lies from equilibrium or else 
produce an equilibrium condition by adding sufficient 
paramagnetic impurity. 

In general, large values of T\ are to be expected for 
S-state atoms because of the weak coupling of the 
electron spin to the rest of the atom and hence, through 
collisions, to the lattice [2]. During the N atom studies 
described in reference [1], H atoms were examined 
briefly [3]. Under some circumstances they were also 
found to be incompletely magnetized. EPR signal 
enhancements of the order 1.3 were observed on addi- 
tion of trace amounts of O2. Recently, this enhancing 
effect on the H atom spectrum has been reported by 
Breen and Glass [4] for trace additions of N0 2 . Westen- 
berg and De Haas [5] have also reported an enhance- 
ment of H atom signals by 2 , but the effect was very 
small and noticeable only at very high carrier flow 



1 Figures in brackets indicate the literature references at the end of this paper. 



rates. Using power saturation methods, they measured 
T\ for H and found values of the order of 100 /xs. These 
are much shorter relaxation times than those observed 
in reference [1] for N atoms. Because their observed 
T\ values were independent of carrier flow velocity 
and increased with increasing carrier pressure, they 
concluded that the predominant relaxation process 
occurred on the flow tube walls. 

The work of Westenberg and De Haas prompted 
the writer to make a more detailed study of spin- 
relaxation effects on the H atom spectrum than that 
reported earlier. Both experimental and theoretical 
aspects of the problem were examined. For the ex- 
periments, H atoms were generated by thermal dis- 
sociation [6] of pure H 2 thereby avoiding the need for 
impurities (such as 2 , N 2 , or H 2 0) [7] normally re- 
quired for efficient dissociation by a discharge. In 
this system traces of 2 added downstream from the 
dissociation point produced significant enhancements 
of the H atom EPR signal. Furthermore, the amount 
of microwave power required to produce a given degree 
of saturation increased with increasing H 2 flow rate. 
These two observations demonstrated that in the 
absence of added 2 , T\ for this system was of the order 
of 10 ms, far larger than that observed in the system 
of Westenberg and De Haas. A simplified analysis 
of the H-H spin-exchange process showed that spin- 
exchange collisions alone can lead to partial mag- 
netization of the atoms in the absence of any added 
paramagnetic impurity. In a typical EPR experiment, 
the maximum H atom signal enhancement to be ex- 
pected on addition of 2 is only a factor of two, in 
contrast to the case of N atoms where much larger 
enhancements are possible. The effects of atom diffu- 
sion on the apparent relaxation time measured by 
power saturation were also examined theoretically. 
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Because of the difficulty in applying corrections for 
diffusion and the additional problems caused by 
inhomogeneously broadened lines and the uncertain- 
ties in estimating the microwave magnetic field 
strength inside the flow tube, no attempt was made to 
measure absolute relaxation times by power saturation 
methods. 

2. Theory 

Hydrogen atoms in their ground electronic state 
possess four magnetic hyperfine levels. In this section 
a. general expression will be derived for the time rate 
of change of the concentrations c, of these levels in a 
cylindrical flow tube. 

A number of processes produce changes in the Cu 
The most important ones are, (a) spin-exchange 
collisions between the H atoms themselves, (b) 
spin-exchange collisions with gas phase paramagnetic 
impurities, (c) collisions with the walls, (d) interactions 
with rf radiation, (e) mass flow, and (f) diffusion. The 
effects of each of these mechanisms will be evaluated 
separately and the results added to give the total rate 
of change in the q. The general equation contains non- 
linear terms arising from the H-H spin-exchange 
process. It will be linearized by considering only small 
departures from equilibrium. To solve the linearized 
equation, further simplifications are necessary. 
Steady state solutions will be obtained by neglecting 
diffusion and assuming a uniform velocity profile 
for the carrier. The effects of axial diffusion will be 
examined by solving a rate equation for a two level 
model, again assuming a uniform carrier velocity. 



2.1. Spin-Exchange Interaction Between H Atoms 

Consider a collision between two H atoms, the 
first in state ^P j and the second in state M^-. We want 
to calculate the probability of the first atom being 
in state ^ and the second in state Wi after the 
collision. In terms of strong field spin eigenfunctions, 
the wave functions ^i, £= 1 — >4, for any value of 
the dc EPR magnetic field, are 






l ' 2 ab 



(1) 



where .a and a are the electron and proton spin func- 
tions respectively, whose angular momentum pro- 
jections in the direction of the dc field have the 
value i, and /3 and b are the functions with projections 
-i; e=[l-x(x 2 +l)-^]l2;x= (I ~ gjlg P ) (vp/vh); 
gj and gp are the electron and proton g factors; vh is 
the zero field hyperfine splitting frequency and 
vp^ — gP/jLoHo/h; Ho is the Bohr magneton and Hq the 
magnitude of the dc field. The energy levels Ui are 
given by the Breit-Rabi equation [8] and lie in the order 
U\ > U 2 > U3 > Ua for a nonzero field. 



Following the treatment of Purcell and Field [9], 
we express the initial spin wave function for the two 
atoms as the product "^(1)^(2). This can be 
written as a superposition of singlet and triplet elec- 
tron spin functions. For example, one might have 

^,7,« ( =^.(l)^2(2)=(l-e)"^(l)a(2)a(l)6(2) 

+ V2€" 2 [a(l)/3(2) + j8(l)a(2)]a(l)o(2) 

+ V 2 e" 2 [a(l)/3(2) -/3(l)a(2)]a(l)a(2). 

The difference in the scattering potential for the 
singlet and triplet components produces a phase 
shift ifts between the amplitudes of the two com- 
ponents. This phase shift is given by the equation 
<fts = j (V t —V s )dt/fi, where V t and V s are the 
energies of the triplet and singlet electronic states, 
and the integration is over the duration of the 
collision. 

After the collision the wave function becomes, 

%inai= (1 - e) " 2 a(l)a(2)a(l)6(2) 

+ 72€" 2 Ml)j8(2) +j8(l)a(2)]a(l)a(2) 

+ 1 / 2 e^[a(l) i 3(2) -j8(l)a(2)]a(l)a(2)exp(^). 

The amplitude of the contribution of a particular 
product function ^(1)^/(2) to ^it final is given by 
the integral / V k (\)W i2)9 finaidr , where the inte- 
gration is over all space. The square of this ampli- 
tude gives the probability that after the collision 
the state ^initial becomes the state ^(1)^(2). 

Performing this calculation for all possible combi- 
nations of initial product functions ^(l) ^(2) and 
all possible final product functions ^(1)^/(2), 
one obtains a 16 X 16 matrix of transition probabilities. 
Because of the small energy differences between the 
hyperfine levels this matrix is not quite symmetric. 
By detailed balancing, the rate of the transition 
ij — » kl, at equilibrium, must equal the rate of the 
reverse process kl — > ij. This can happen only if the 
matrix elements satisfy the relation (kl, ij) = 
(ij, kl) exp [(U k + Ui-Ui-U 3 )lkT], where (ij, kl) 
is the probability for the transition ij — > kl, k is Boltz- 
mann's constant, and Tis the gas kinetic temperature. 
The error involved in calculating the transition prob- 
abilities in this manner is discussed by Purcell and 
Field. It is of the order (U k + Ui-Ui-Uj)lkT and 
thus quite small. 

The nonzero elements of this matrix are given in 
table 1 as a function of the field strength parameter 
e and a phase shift parameter A = i sin 2 (<p^/2). Only 
strong collisions for which <pt s > 1 will be considered. 
In this case A = i. Transition probabilities for the zero 
field case treated by Purcell and Field can be obtained 
from table 1 by letting e = i. For an EPR absorption at 
X band frequencies, e ~ 0.5 X 10~ 2 for the high field 
H atom line. The matrix element (ij, kl) is only the 
relative probability that y— > kl in a spin-exchange 
collision. The overall rate at which such collisions 
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occur will be specified by a rate constant S. A litera- 
ture value will be used for this parameter. 
The rates of change of the c* are given by 

c<S ~* = - 2 aq - 2 Ciq 



+ 2 c ^ { 2 («. w») + 2 («. ™) } 



A:,/ 



m#/ 



: ^ a k nc k ci, 
k, I 



;=i- 



(2) 



where the coefficients a^ are various linear combina- 
tions of the matrix elements (i/, A:/). These equations 
represent changes in the a arising from spin-exchange 
collisions between the H atoms themselves. Let us 
now consider changes resulting from interactions of 
the atoms with the external system. 

TABLE 1. Nonzero transition probabilities for spin-exchange 
collisions between H atoms, (ij, kl). 



Transition ■ 



(12, 14), (12,41) 

(13,22), (13,44) 

(14,21) 

(21,41) 

(22, 31) 

(23, 34), (23,43) 

(31, 44) 

(32, 34), (32,43) 

(12, 21) 
(13, 24) 
(34, 43) 

(13,42) 
(14,41) 
(23, 32) 

(22, 24), (22,42) 
(24, 44) 
(42, 44) 

(22, 44) 
(24, 42) 

(11, 11), (33,33) 

(12, 12), (21,21) 
(34,34), (43,43) 

(13, 13), (31,31) 

(14, 14), (41,41) 
(32, 32), (23,23) 

(22, 22), (44,44) 



(24, 24), (42,42) 



Transition probability 



e(l-e)A 

e-A 
(1-€)*A 

6(1-6) (1-2€)*A 

46 2 (1-6) 2 A 

1 
l-2€(l-i€)A 

1-A 

1-(1-€*)A 



1-46(1-6) 

[1-6(1"6)]A 

l-[l-46-(l-6) 2 ]A 



{kl,ij)=(ij\kl)exp [(U k + Ui-U t -Uj)lkT] 



a The matrix element ((/, kl) corresponds to the transition 
i//,(l)i|/j(2) -> i//a-(1)i///(2), where the numbers in parentheses label 
the two atoms involved in the collision. 



2.2. Spin-Lattice Interactions 

For H atoms in a nonpolar, nonparamagnetic carrier, 
the only important spin-lattice relaxation processes 
are spin-exchange collisions with paramagnetic impuri- 
ties and wall collisions. The paramagnetic impurity is 
assumed to have orbital angular momentum so that its 
electron spin will be strongly coupled to the lattice. The 
spin temperature of this species will thus always 
remain at the temperature of the lattice [2]. We shall 
not try to specify the details of the impurity interactions, 
but shall simply let Kq be an overall rate constant for a 
particular transition Wr+tyj produced by these 
interactions. 

The total rate of change in the c* resulting from these 
transitions will be 



c\ = ^ (KkiCk — Kikd) = ]T RijCj 



(3) 



where Rij = Kji and /?,•/ = — \ Kij. From detailed 

balancing we have Kij^Kjj exp [(U i — V j) I kT\ One 
would expect the Kq to be proportional to the con- 
centration of a particular paramagnetic impurity and to 
depend on its nature. 

Relaxation by wall collisions also yields terms which 
are first order in the c,. If these were treated together 
with the gas phase interactions a single set of Kq 
parameters could be used to represent both processes. 
It will be shown later that the wall relaxation process 
can be included explicitly by treating the interaction as 
a boundary condition on the rate equation when radial 
diffusion of the atoms is included. 

2.3. Transitions Induced by an rf Magnetic Field 

The rate of change in the ct arising from absorption 
and induced emission by an rf magnetic field is given 



k*i j 



(4) 



where the Wq are rate constants for these transitions, 
w\j=Wju and wn = — ^Wij. In the present experi- 

ment, the high field H atom line which arises from tran- 
sitions between levels 2 and 3 was examined. For 
this case all Wij = except W23 and Wyi which are 
equal. Let W 2 z and W32 equal W. The rate constant W 
is proportional to the rf power striking the sample. Be- 
cause the rf field varies from point to point throughout 
the microwave cavity, W will be different at different 
points 'in the cavity. We have W=h-' 1 H\fx 1 g{v), 
where H x is the amplitude of the rf magnetic field, fi 2 
is the square of the absolute value of the transition 
matrix element, and g{v) is the normalized line shape 
function [10]. 

2.4. Effects of Mass Flow and Diffusion 

To complete the rate equations for C|, one must 
include a term which accounts for the flow of atoms 
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through the microwave cavity and another which rep- 
resents their diffusion. The experimental arrange- 
ment consisted of a cylindrical flow tube running 
through a TEo n cavity with its axis coinciding with 
that of the cavity. Were relaxation on the walls as 
effective a 1 in 10 collisions in the present experi- 
ments, large radial gradients would arise in the pop- 
ulation differences between the levels of interest 
during a saturation experiment. In this case both radial 
and axial diffusion terms should be included. 

2.5. Complete Rate Equations and Their Linearization 

The complete rate equations for the c, are, in cylin- 
drical coordinates 



c)Cj 



-TT=S 2 a kU c k ci + ^ RijCj + ^j Wi > c 



dt 



k, I 



- v{ r)¥ + D\±±(rf)+%}} (5) 

dz [ r ()r\ or ) dz 1 J 

where v{r) =2vo[l— (r/a) 2 ] the carrier flow velocity, 
vo is the average carrier velocity, a is the radius of the 
flow tube, and D is the atom's diffusion coefficient at 
the prevailing carrier pressure. The first term in the 
brackets represents radial diffusion and the second 
axial diffusion. The term containing v(r) accounts for 
the effects of mass flow. For reasons discussed earlier, 
the coefficients Wij, which represent the effect of rf 
power, are functions of r and z. The sum containing 
the parameters Rij represents the effect of interactions 
with paramagnetic impurities and if desired, the effect 
of relaxation on the walls. The first sum represents the 
effect of spin-exchange collisions among the H atoms 
themselves. 

Under all conditions experienced in a typical EPR 
experiment with H atoms, the c, differ at most by 
several parts per thousand. The rate equations there- 
fore can be linearized [11] by expressing the c\ in the 
form c/=ira(H-fi), where n is the total H atom con- 
centration and the f ,• are very small numbers which 
will be considered only to first order. As a result of 
this definition, the relation V £/ = must always be 

i 

satisfied. 

With this substitution, the rate equations become. 



dt 



k, I 



y JbJ v J dz [r dr\ dr J dz 2 J 

where the relation V Wij = has been used. The Ay 

j 
are linear combinations of the a^u, which in turn are 
linear combinations of the transition probabilities 
(ij, kl) from table 1. 

The Boltzmann factors in the (ij, kl) and in the rate 
constants K\j from which the Rij are derived ensure 



that these rate equations will satisfy the proper 
equilibrium conditions. However, the equations are 
easier to handle if they are rewritten so that the 
equilibrium values of the £; are included explicitly. 
This may be done in the following manner. For the 
system at thermal equilibrium, with no rf power 
applied, the rate equations become 



I 



1 



Sn X a kli + X *«=-7 Sn J A^j-^R^j 



where the E\ are the equilibrium values of t;\ for the 
particular lattice temperature and dc field. 
Substituting this expression into (6) gives 



d| i= l 

dt 4 



Sn^Auitj-F^+^Rui^-Ej) 



+Wj -„ w f +0 {I|( r f) + SM 



In this form the equations give the correct equilibrium 
values of £, even if the Boltzmann factors inherent in 
the Aij and Rij are neglected. The error in doing so 
will be very small since the factors are so close to 
unity. 

2.6. Evaluation of the Coefficients A- t j and Rq 

The coefficients A\j which account for the effects 
of H-H spin-exchange collisions are linear combina- 
tions of the transition probabilities (ij, kl) shown in 
table 1. They are given in table 2 in the form of a 
4X4 matrix which is valid for all values of the dc 
field. Strong collisions, for which A = i, have been 
assumed, and the Boltzmann factors neglected. In 
this table an additional field strength parameter, 
X = 8e(l — e) has been introduced. Lambe [11] has 
also determined the A\j coefficients by a somewhat 
different method. 

In the strong field limit, the coefficients Rij which 
describe the spin-lattice interactions can easily be 
derived if we consider only interactions which change 
the electron spin state. These interactions could be 
with a gaseous paramagnetic impurity or with the 
walls of the flow tube. On changing the electron spin 
quantum number, the strong field wave functions for 



Table 2. The coefficients A u which describe the effects of H-H spin- 
exchange collisions. 

The parameter x equals 8e(l — e). 



i J 


1 


2 


3 


4 


1 


-1 


1 


-1 


1 


2 


1 


-d-r-X) 


1 


-d-x) 


3 


-1 


1 


-1 


1 


4 


1 


-d-x) 


1 


-d+x) 
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the four hyperfine states transform in the following 
way: 

^ 2 = ab -> fib = W :] 

The rates of change in the c, produced by this 
interaction are thus 

C\ = K41C4 —K14C1 
C2 = K32C3 — K23C2 

C:i = K 2 3C2 — K32C3 

c a = K\aC\ —K 4 \Ca. 



where Kjj = iK, and K is an overall rate constant for 
these interactions. The reason why the factor i has 
been introduced will be discussed later. The Rij can 
be derived from the K\j by means of (3). 

By considering only the strong field limit, one is 
assuming that the spin-lattice interaction effects only 
the electron spin. Although this leads to equilibration of 
the population differences C\ — c 4 and Co — c-.\ which are 
the ones involved in the EPR transitions, it will not 
by itself produce thermal equilibration of the d. This 
is because transitions between the states ^1, ^4 and 
^2, ty-A have been neglected. Such transitions involve 
changes in the proton spin state. One way to estimate 
the relative Kq for such changes is through an adapta- 
tion of the method used by Lambe [11 J to determine the 
rate equations for H-H spin-exchange collisions. 
Consider the paramagnetic impurities to be a reservoir 
of electron spins in which P* and P~ are the probabili- 
ties of finding an electron with spin up and spin down, 
respectively. Since the impurity spins are assumed to 
be strongly coupled to the lattice, and thus always at 
thermal equilibrium, P f and P~ will be constant and 
have almost the same value, namely i, at the temper- 
atures and field strengths used in the present experi- 
ments. Suppose that an H atom in state ^1 exchanges 
its electron with one from the impurity reservoir. 
With probability / , + , ^1 remains unchanged. With 
probability P~, ^1 goes to state /3a which is a mixed 
state. We have 

/3a= (l-eV^ + e 1 ^. 

so that with a probability eP~, ^1 goes to "^2, and with a 
probability ( 1 - e)P , to %. 

Following this procedure for all the states ^/, one 
gets rate equations for the c\. For example, 

c ! = -KPcx + KePc> + £( 1 - e)P + c 4 
= -i/<c, +iKec 2 +i#(l-e)c 4 

where K is an overall rate constant for spin-exchange 
collisions between the H atoms and the paramagnetic 



impurity. It will be proportional to the impurity 
concentration. 

On making the transformation c, =in(\ + f/), and 
introducing the equilibrium parameters Eu one ob- 
tains the rate equations in the desired form f, 
= ^jRij{% :) —Ej), where the matrix Rij is given in 

J 
table 3. For the strong field case, e— »0, and the Rij 
reduce to the values obtained earlier. The factor \ 
in these coefficients arises because \ of the impurity 
spins in the strong field case will be in the same direc- 
tion as the electron spin on the H atom and so no 
change in the state of H will occur during an electron 
exchange in one out of two collisions. 

While the matrix R l3 in table 3 should be valid for 
all field strengths, only the strong field elements were 
used to analyze the experimental data. The H-H spin- 
exchange interactions also produce the transitions ne- 
glected in the strong field spin-lattice approximation. 
Since the H-H exchange dominates the spin-lattice 
interactions in the present experiments, the spin- 
lattice contributions to these transitions can be ignored. 

Table 3. The coefficients R/j, which describe the spin-lattice 
interactions. 

The parameter x equals 8e(l — e). 



i \ 


1 


2 


3 


4 


1 

2 
3 
4 


-iK 

UK 



kK 
i(l-€)K 

xW 




i(l-€)A 

-iK 

UK 


h(l-e)K 
XK/S 
kK 

-i(i+ix)£ 



2.7. Wall Relaxation as a Boundary Condition 

If relaxation on the walls were rapid enough to 
produce radial concentration gradients, the rate 
equation (7) should be solved subject to the boundary 
conditions 



$1 -4«!> 2 (*-M6-*i> 



where the derivatives are evaluated at the wall of the 
flow tube, D is the diffusion coefficient for the atom in 
the carrier gas at the prevailing carrier pressure, and 
the (Rw)ij are coefficients for wall relaxation analogous 
to the Rjj defined for relaxation by paramagnetic im- 
purities in the gas phase. In the present experiments 
wall relaxation cannot have been rapid. Therefore the 
radial diffusion term in (7) will be dropped, and any 
wall relaxation effects will be absorbed in the Rij 
coefficients instead of being included explicitly as a 
boundary condition. 

2.8. A Two-Level Spin Model 

In order to examine the effects of axial diffusion, a 
two-level model consisting of particles with spin V2 was 
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used. For such a system, the Aij = 0, while the matrices 
for Rij and wij are simply 



Rtj z 



IK 






and, 



Wu= 



W W 
W -W 



(8a) 



(8b) 



where K and W are analogous to the same quantities 
defined for the four-level model. 

3. Apparatus and Experiments 

The apparatus used was essentially that described 
in reference 1 except that the klystron frequency was 
now locked to the cavity resonance instead of to an 
external oscillator. Field modulation was at 400 Hz. 
Hydrogen atoms were produced outside the high 
field of the EPR magnet by passing Matheson 2 Ultra 
High Purity H2 over a tungsten filament [6] heated 
to about 1600 °C. Only H2 was used as a carrier. The 
walls of the 10 mm in diam quartz flow tube were 
poisoned against atom recombination with boric 
acid [12]. The amount of atom decay was measured by 
adding a trace of NO to the atom stream and observing 
the HNO emission [13]. Between the EPR cavity and 
the point of O2 addition, the decay was less than 
5 percent for all experiments. The length of flow tube 
in the high-field zone upstream from the cavity was 
estimated to be about 15 cm. Impurity O2 was normally 
added at a point 45 cm downstream from the filament 
and 27 cm upstream from the cavity center in the 
form of a 1:500, 2 : Ar mixture. At H 2 flow rates of 
the order 3.5 atm cm 3 s _1 , and pressures of 2 torr, a 
2-fold enhancement of the unsaturated high field H 
atom EPR line was observed on addition of 15 ppm of 
O2. Oxygen additions greater than this produced no 
further increase in the signal, nor did the addition of 
this same 15 ppm of O2 further upstream. 

At constant atom concentrations in the cavity, the 
amount of power required to produce a given degree 
of saturation was found to increase with increasing 
carrier flow velocity. Carrier velocities from 650 to 
4400 cm s _1 were covered in these saturation 
experiments. 

Analysis of these experiments indicated that 7\ 
in our system was of the order of 20 ms in the absence 
of added O2, a value similar to that observed for 
N atoms. The neglect of the radial diffusion term in 
(7) is therefore justified, since even if a wall mechanism 
were the sole relaxation process, the radial gradients 
it would produce would be negligible at the flows and 
pressures used. 



2 Certain commercial products and instruments are identified in this paper in order to 
specify adequately the experimental procedure. In no case does such identification imply 
recommendation or endorsement by the National Bureau of Standards, nor does it imply 
that the products or equipment identified are necessarily the best available for the purpose. 



(9) 



4. Solution of the Rate Equations 

Since a flow method was used, only steady-state 
solutions of the rate equations were required. The 
equations (7) for the four-level system were solved 
only under the following conditions: (a) radial and axial 
diffusion neglected; (b) for the unsaturated case, a 
uniform velocity profile for the carrier, and (c) for the 
saturated case, a zero carrier velocity. To examine the 
effects of axial diffusion at low and high power, and 
nonzero carrier velocity under saturation conditions, 
a two-level model was analyzed. Again radial diffusion 
was neglected, and a uniform carrier velocity assumed. 

4.1. Unsaturated Four-Level System 

For low power conditions such that the coefficients 
wij are negligible, diffusion is neglected, and the 
carrier velocity is assumed to be uniform, the steady- 
state solution of (7) is given by the set of equations, 

fr = E 3 +i(V-Y) 

£ 4 =E 4 +i(U+V) 

where, 

F=F„exp [-2(2y+lK) (zM] 

Y=H- l {(Y u e + -iKU ) exp(<D-) 

+ (Y,e- + iKU„)exp(^)} 

U=U- % {(U e--iKYo) exp(<p-) 

+ (U d + +iKY„)exp (<*>+)} 
and, 

y — iSn 

t=(r)*-Kyx) 112 

± =i±7X 

<t>± = -( V ±0(z/vo). 

The subscript zero on the parameters U, V, and Y 
denotes their values at z=0, the start of the high field 
region. As mentioned earlier, the strong field values of 
the Rij coefficients were used. 

Calculated values of the f * parameters are shown in 
figures 1 and 2 as functions of z/vq. For this calculation 
the gi values at the start of the high field region were 
taken to be the values corresponding to thermal 
equilibrium at zero field and 300 K, namely, £? = f° 
= £> = - 0.572 XlO" 4 , and fS= 1.716 X10- 4 . An H atom 
concentration of 0.44X10 14 cm -3 and a H — H spin- 
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FIGURE 1. Population parameters £, versus z/vo calculated 
from eq 9. 

Parameter values used were K = 0, S=13X10~'° cm 3 s~', n = 0.44XlO M cm -3 , dc 
field = 3500 Gauss, and T=300 K. The difference &-£ 3 is also plotted. 
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Figure 2. Population parameters £, versus z/v calculated 
from eq 9. 

Parameter values used were K= 168 s~\ S=13xl0-'° s~\ n = 0.44Xl0 14 cm" 3 , dc 
field - 3500 Gauss, and T= 300 K. The difference £, - £ 3 is also plotted. 



exchange rate constant of 13 X 10 -lo cm 3 s -1 were used. 
In the present apparatus the EPR cavity would lie at 
z/vo= 12XlO- 3 s for *; = 800 cm s" 1 . In figure 1, the 
spin-lattice relaxation parameter K has been set equal 
to zero, while in figure 2, it was given a value of 168 s _1 . 
Also shown in both figures is the difference £2 — f.3- 
The intensity of the high field EPR line will be propor- 
tional to the value of this difference at the cavity posi- 
tion. Thermal equilibrium values of the parameters at 
a field of 3500 Gauss and 300 K are shown as dotted 
lines. 

In the absence of any spin-lattice relaxation, £1 and 
£3 remain unchanged, while H-H spin-exchange brings 
about rapid relaxation of £ 2 and £4 to their equilibrium 
values - 8.431 X 10~ 4 and 8.490 X 10~ 4 , respectively. 
They approach to within 10 percent of their final 
values in a time of the order l/3yx- Because x 1S so 
small (~ 0.04), this is considerably slower than one 
might expect from the large value of y at this particular 
H atom concentration. It is obvious that the strong 
field wave functions would not be suitable for this prob- 
lem because they yield x~0 an d tnus no relaxation of 
£2 and £ 4 . 

The absolute value of £ 2 — £3 for this partially re- 
laxed condition is 6.77X10 -4 , which is a factor of 2.16 
smaller than the thermal equilibrium value. Therefore, 
unless one operates at considerably smaller H atom 
concentrations, or at much higher flow rates, the maxi- 
mum enhancement of the EPR signal to be expected 
on addition of a paramagnetic impurity would be this 
factor of 2. 16. 



In case y\ > K, the solution for ^2 — £3 may be 
written approximately as 



£2-e3=E 2 -Es+i{Y -Uoexv(-2yxzlvo)} 

exp(- 



iKz/vo). 



The term in brackets quickly relaxes to the value Yo 
so that the process can be described by a single relaxa- 
tion rate over the range of slow decay. 

Let us now assume that K may be written as a sum 
of two terms, i.e., £=£o + (?[02], where Q is a rate 
constant for spin-exchange collisions between H and 
O2, [O2] is the O2 concentration, and Ko is the value 
of K without added O2, presumably arising from colli- 
sions with paramagnetic impurities and/or wall relaxa- 
tion. From the observed values of the high-field line 
intensity as a function of added O2, and the results of 
calculated £, like those shown in figures 1 and 2, a 
value of 6.4±2.0 X 10~ 10 cm 3 s _1 (corrected for axial 
diffusion as described in the next section) was found 
for Q, and a value of 50±20 s -1 for K at a H 2 flow of 
3.5 atm cm :i s _1 and a pressure of 2.2 torr. No measure- 
ments were made at other flows or pressures. 



4.2. Unsaturated Two-Level System and Corrections 
for Axial Diffusion 

For a two level system, a steady-state analytical 
solution of (7) is easily obtained if radial diffusion is 
neglected and a uniform carrier velocity is assumed. 
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The solution for the population difference 8£ between where 
the two levels in the unsaturated case is 



8f = 8£+(8f°-8£) exp \ -^-rjz —z 



where 8E is 8£ at equilibrium and 8£° is 8£ at 2 = 0, the 
start of the high field region; 8£° will be assumed to be 
zero. 

If vl > iKD, the solution may be approximated as 

S£~ 8E+8E exp {- K[l - (KD/vi)] (z/vo) }. 

This result was used to estimate the effects of axial 
diffusion in the four-level system. For the O2 signal- 
enhancement experiment discussed in the preceding 
section, the correction for this axial diffusion to the 
rate constant Q was small, approximately + 10 percent. 

4.3. Four-Level System With Saturation 

Consider equations (7) at high microwave power 
levels where the factors w\j cannot be neglected. For 
a nonflowing gas, without diffusion, the steady-state 
solution for the population difference of interest is 

h-b^iEt-EzHl + TtW)- 1 

where T\ is given by the expression 

T,=K^{Kly + 2){\Kly + 2)-K (10) 

This parameter is equivalent to the reciprocal of the 
relaxation probability parameter W r defined by Lloyd 
and Pake [10]. The factor multiplying K~ x ranges 
between the limits 1 and 2 according to whether the 
ratio K/y is small or large. Although K/y is much less 
than unity in the present experiments, large O2 addi- 
tions can easily make K > y. In this situation, T\ will 
not be strictly proportional to [O2] 1 . It also is evident 
that when /C/y~l, 7\ will be somewhat dependent, 
through 7, on the H atom concentration. 

4.4. Two-Level System With Saturation 

To estimate the effect of flow and axial diffusion on 
the amount of microwave power required to produce 
a given degree of saturation, the two-level model has 
been used. The parameter W has been assumed to 
be constant over the effective length L of the cavity 
(~ 2cm) and zero elsewhere. The spin-lattice parame- 
ter K has been assumed to be large enough to equil- 
ibrate the spins by the time the atoms reach the 
cavity. The EPR signal will be proportional to the aver- 
age value of the population difference 8£ in the cavity. 
This average over the length of the cavity is given by 



F± = 2(C ± //x ± L) sinh (fiHL) 

C--=±WdE(iK+W)-*{(^- f jL+) f jL-exv(- f jLHL) 

-(/i--ju,o)/4 exp (/xHI)}(G + +G-)- 1 
G ± =±(/x + -^)(/x--^ ( |)exp [±(^ -fi-)iL] 
/x ± = vol2D ± X 

fJL^ = Vol2D±Xo 

\=[v* + 8D(iK+W)yi 2 l2D 

\o = [v* + 8DiKyi*l2D. 

This rather formidable expression was used to 
calculate (8£) as a function of W for various values of 
the parameters, vq, D, and K. It was found that the 
relationship between (8£) and W could be represented 
to a very good approximation by the expression 



<8£> = s£(i + r*r)- 



(12) 



(8£) = F + + F-+iK8E(iK+W)- 



(11) 



where T* depends on u , D, and K. 

The relaxation time for the four-level model at 
high H atom concentrations will actually be equal to 
jrf. This is because H-H spin-exchange collisions 
couple in the two unpumped levels which are also 
experiencing the same spin-lattice relaxation processes. 

Some examples of how iT* depends on the various 
parameters are shown in figures 3 and 4. Figure 3 
shows the variation of 2/T* with D for a range of K 
values from 100 to 1600 s _1 . The flow velocity vo 
was varied over the range 100 to 500 cm s _1 Over 
this interval the dependence on vo was small and is 
indicated by the width of the curves. It is evident that 
the effect of diffusion can be quite large even in the 
absence of any wall relaxation. 

The variation of 2/7 1 ,* with v Q at ,D = 912 cm 2 s" 1 is 
shown is figure 4 for different values of K. It can be 
seen that for K values above 400 s _1 , the dependence 
of 2/r* on vo is quite small for vo in the range to 1000 
cm s _1 . 



4.5. Determination of 1\ by Saturation Methods 

The fact that the population differences of the 
pumped levels for the models considered here are 
proportional to a quantity of the form {\ + T\W)~ x 
means that the saturation theory for Lorentzian lines 
discussed by Westenberg and De Haas [5 J is applicable. 
Caution must be observed in using this theory, how- 
ever, as it applies only to homogeneously broadened 
lines and for field modulation frequencies co m such 
that cx) m ^ T~ l . The 400 Hz modulation frequency 
used in the present experiment satisfied the latter 
condition when the carrier flow was sufficiently high, 
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D (I0 3 cm 2 sec 1 ) 



Figure 3. Variation of 2/T : 



D calculated from cijs II and 



12 for different values ofK. 

For each K value, v was varied over the range 100 to S00 c 
2/7* is indicated hy the width of each curve. 



is '; the effect of this on 



but the former requirement was not satisfied since 
the lines were asymmetrical and presumably broad- 
ened by field inhomogeneities. 

4.6. Results of Saturation Experiment 

The signal enhancement produced by traces of 
2 demonstrated that the H atom spin-lattice relaxa- 
tion time in the present system was quite long. In this 
case, as can be seen from figure 4, the relaxation time 
observed in a saturation experiment should depend 
markedly on the rate at which unsaturated atoms enter 
the cavity. The effect of carrier velocity on the amount 
of saturation produced by a given microwave power 
level is shown in figure 5. Here, a H 2 carrier was used 
at 2 torr with v = 1160 and 3239 cm s _1 . The high fre- 
quency H atom line is shown at a low power level at 
which saturation was negligible, and at a power level 
20 dB higher for which considerable saturation was 
evident. In the double attenuator method used [1], the 
signal is proportional to x'\ the imaginary part of the 
magnetic susceptibility. It is obvious that there was 
considerably less saturation at the higher power level 
for the higher carrier flow. In this experiment, about 
40 ppm O2 was added to the carrier to ensure complete 
relaxation of the atoms arriving at the cavity so that 
the signal observed at low power would be propor- 
tional to the H atom concentration. This concentration 
was the same for both carrier velocities. The con- 
ditions of this experiment correspond to the third 




12 3 4 

v o (l0 3 cm sec-') 

FIGURE 4. Variation of '2/Tf with v at D = 912 cm 2 s~ l , calculated 
from eqs 11 and 12 for different values ofK. 



v= 1160 cm 



v=3239 




FIGURE 5. Effect of carrier flow velocity on the power saturation of 
the high-field H atom EPR line. 

The two traces labeJed (a) were recorded at low power where saturation was negligible. 
The H atom concentration is the same in both traces but the carrier velocity was different 
as shown. The traces (b) were recorded at a power level 20 dB higher. The lower trace in 
(b) corresponding to the higher carrier flow is less saturated than the one at the lower flow. 
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curve in figure 4, labeled £—1000 s" 1 . This calcula- 
tion predicts that a given power input would produce 
roughly half the saturation at i>o = 3239 cm s _1 that it 
would produce at 1160 cm s _1 . The results shown in 
figure 5 are in agreement with this prediction. 

5. Discussion 

The purpose of this work has been to show how the H 
atom hyperfine populations equilibrate when the atom 
is placed in a magnetic field. In the present experi- 
ments there were two major relaxation processes, (a) 
spin-exchange collisions between the H atoms, and (b) 
spin-exchange collisions between the H atoms and 2 
molecules added in trace amounts. The spin-exchange 
process conserves the component of angular momentum 
in the direction of the dc magnetic field, which we 
shall call the z-direction. Because states ^ 2 and ^ 4 
have z angular momentum components of zero, spin- 
exchange collisions among the H atoms themselves 
can equilibrate [9] these levels without changing the 
net value of the z-component. The net magnetization of 
the system will change however, since these two 
states have nonzero magnetic moments at high field 
strengths. Thus, partial magnetization of the system 
will arise solely from H-H collisions. It is worth noting 
that this situation cannot occur for N atoms, because 
none of the hyperfine levels have zero values for the 
z-component. 

The other two states have nonzero components 
-\-h and — h for ^i and ^ 3 , respectively. Thus, the 
relative populations of these two levels cannot change 
without changing the total value of the z-component 
for the H atom system. This can be changed only 
through an interaction with the surrounding system, 
i.e., with the "lattice". Spin-exchange collisions with 
impurity 2 can provide the necessary spin-lattice 
interaction even though the z-component will be con- 
served during the collision. This is because any non- 
thermal spin orientation imparted to the 2 will be 
rapidly destroyed because of the strong coupling of the 
electron spin on 2 to its molecular framework. This 
is strongly coupled to the heat bath particles through 
collisions which change the 2 rotational energy. 

From the enhancing effect of 2 on the H atom 
signal, a rate constant of 6 X 10" 10 cm 3 s" 1 was deter- 
mined for the H-0 2 spin-exchange collision. Kleppner, 
Goldenberg, and Ramsey [14] examined the effect of 
2 on the radiating state of the atoms in a hydrogen 
maser, and reported a relaxation rate constant. To 
compare our value with theirs, it is necessary to solve 
(7) for the zero-field case. Letting e=i in tables 1 and 
2 for Aij and Rij, and neglecting diffusion and mass 
flow, one obtains for the population difference f 2 _ ^4 
involved in the maser transition, 



6-£4=£ 2 -£4+(£ 2 °- 



-£ 4 °-£ 2 +£ 4 )exp 



will be represented by the change of exp (— Kt) with 
time. Therefore the relaxation rate observed by these 
workers is equal to K as defined in the present work. 
Their value of K[0 2 ]- l = 6.2 X lO" 10 cm 3 s" 1 is in 
agreement with ours. Neither of these values is in 
agreement with that of Westenberg and De Haas [5], 
whose t\ values should be equal to K~ l since their 
H atom concentration was kept high during 2 addition. 
(See eq. 10.) 

A collision cross section cr 2 can be derived from the 
observed rate parameter K^Q[Oz] if one uses the 
kinetic-theory relation K=cr' 2 (8irkT/iJL) 1/2 [0 2 ], where 
k is the gas constant, T is the absolute temperature, 
and jx is the reduced mass of the colliding particles. 
The values of the rate constant Q and the cross section 
cr- are shown in table 4. 



Table 4. Rate constant and collision cross section for H— O2 
spin-exchange collisions. 


Q 


& 1 


Reference 


6.4±2.0x 10- ,o cm :{ s- 1 
6.2 
35 


8.0±2.5X10- I,: cm 2 

7.7 
44 


This work. 

[12] 

[5] 



(- 4,y t-Kt) 



where t is the time, and f 9 are the values of & at t = 0. 
If 4>y< K in the presence of added 2 , the relaxation 



The cross section determined in the present work 
and by Kleppner, et al., is smaller than that for H-H 
spin-exchange collisions. This may not be unreasonable 
since one would expect the attractive potential curve 
for H~0 2 to be Considerably shallower than that for 
H-H and thus the phase shift between the attractive 
and repulsive potentials could be smaller. However, 
since it appears that the steepness of the repulsive 
state potential curve is also an important factor, at 
least in the H-H case [15], it is difficult to say anything 
definite about the H-0 2 system without a detailed 
knowledge of the states involved. 

It is not clear why wall relaxation should have been 
so much more important in the experiments of Westen- 
berg and De Haas than in ours. Nevertheless, our 
observations show that it is quite possible for wall 
relaxation to be negligible in some systems. In this 
situation some 2 impurity must be present to achieve 
complete magnetization. The amount necessary at 
moderate carrier flows is very small however, of the 
order 15 to 30 ppm, so that it should not interfere with 
kinetic studies of moderately fast reactions. Under 
such conditions however, power saturation can still 
depend on the carrier flow rate. A line which is un- 
saturated at a high flow rate may be saturated at a 
lower one. 
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